NASA Contractor Report 195421 


//!/- 7/ 

w 


The Use of Cowl Camber and Taper to Reduce 
Rotor/Stator Interaction Noise 


R. Martinez 

Cambridge Acoustical Associates, Inc. 
Medford, Massachusetts 


(NASA— CR-195421) THE USE OF COWL N95-22675 

CAMBER AND TAPER TO REDUCE 
ROTOR/STATOR INTERACTION NGISE 

Final Report (Cambridge Acoustical Unclas 

Associates) 89 p 


G3/71 0042157 


February 1995 


Prepared for 

Lewis Research Center 

Under Contract NAS 3-27229 



National Aeronautics and 
Space Administration 



CAMBRIDGE ACOUSTICAL ASSOCIATES, INC. 


THE USE OF COWL CAMBER AND TAPER TO REDUCE ROTOR/STATOR 

INTERACTION NOISE 


by R. Martinez 


TABLE OF CONTENTS 


Page 


I ABSTRACT 

H INTRODUCTION * 

m INCIDENT FIELD ANALYSIS 

A Chordwise-Compact and Chordwise Noncompact Source Representations for 
Twisted Stator and Rotor Blades * 

1 . General 

2. Description of a Three-Dimensional Stator/Rotor Geometry for the 

Generation of the Cowl's Internal Incident Neaifield 

3. Mathematical Construction of the Impinging Wakes 

4. The Radiated Nearfield for both Chordwise-Compact and Noncompact 

Blade Loadings ; 

5. Chordwise-Compact and Noncompact Blade Airload Solutions 


B Farfield Expressions for p * ^ 

1. General 

2. Farfield for Rotor Sources in the Aircraft's Reference Frame 

3. Incident Farfield for the Ground Reference Frame - 

4. Conversion of the Blade-Source Analysis in Sections A.2-B.3 to Apply to a 

Rotor-First Stage 


IV A CANONICAL MODEL OF COMBINED LEADING- AND TRAILING-EDGE 
DIFFRACTION FOR THE COWL: 

ACOUSTIC SHIELDING BY AN UNWRAPPED DUCT 

/ 

A Motivation ‘ 

B Brief Background of Screen Diffraction Problems with and without Flow ! 



CAMBRIDGE ACOUSTICAL ASSOCIATES, INC. 


C Postulated Reconciliation of Existing Leading-and Trailing-Edge Canonical 
Solutions ••••’• 25 

D High-Frequency Justification of a Combined Leading- and Trailing-Edge 

Diffraction Model 27 

E Numerical Predictions of Noise Shielding by the Canonical Model: Effects of 
Cowl Camber and Taper 27 

V SAMPLE PREDICTIONS OF FORMAL COWL DIFFRACTION FOR 
STATOR-FIRST, AND ROTOR-FIRST SYSTEMS 

A Kinematics of Ray tilting in a Uniform Freestream 30 

B Rigorous Predictions of Cowl Diffraction for Section IILB's Incident Fields .. 31 

1. The Unlined Cowl — 31 

2. The Lined Cowl 35 

YI CONCLUSIONS 38 

ACKNOWLEDGMENTS 38 

REFERENCES 39 

FIGURES 41-66 

APPENDIX A: COUPLED LIFTING- AND NONLIFTING- SU RFACE THEORIES 

FOR A STRAIGHT COWL A1 ' 6 

APPENDIX B: COUPLED LIFTING- AND NONLIFTING-SURFACE THEORIES 

FOR A COWL OF ARBITRARY CAMBER AND TAPER B. 1 -9 

PROJECT DOCUMENTATION PAGE - c 



CAMBRIDGE ACOUSTICAL ASSOCIATES, INC. 


NOMENCLATURE 


Roman 

a cowl radius normalized by the duct's halflength L/2 

A parameter used to characterize the width of the 

impinging wakes 

B r number of rotor blades 

B s number of stator blades 

c speed of sound of the medium 

c(r *) nondimensional chord for the radiating blades 

D/Dt substantial derivative, dldt + U m d/dz 

H, Hankel function of the first kind of order v 

J v Bessel function of order v 

k cowl's reduced frequency, coL^U^ 

k transform wavenumber; nondimensional in Appendix A and dimensional in 

Appendix B 

L/2 cowl's halflength or halfchord, the normalizing constant for all distances and 

wavenumbers 


L Uner length of lined segment (Appendix A) 

m temporal harmonic counter; m=l corresponds to the blade-passage frequency; 

also serves in the same context as the rotor wake’s Fourier counter 

M ax flight Mach number; the "ax" subscript denotes "axial" 
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Mach number relative to the blade section at radial station r* 

Stator Fourier counter 

time-dependent pressure; dimensional; as on the left-hand side of Eq. 5a 

frequency-domain pressure, and its circumferential v mode; both dimensional 

incident field pressure of circumferential mode v , corresponding to exp (-zcof) . 

lifting and thickness contributions to the scattered pressure; dimensional 

pressure source strength in the traditional treatment of the trailing-edge diffraction 
problem 

volume source strength in the canonical model 

field value of the radial coordinate; r=a corresponds to the cowl's wall; 
nondimensional 

running radial variable distinguishing among blade stations; nondimensional 

nondimensional range to the farfield observer, with origin at the cowl's geometric 
center r=0, z=0. R q is its dimensional version, RoL/2, used in the retarded time 
on Eq. 16's left-hand side. 

distance from vertex to source point in the canonical model; dimensional 
Sears function 
time; dimensional 

aerodynamic transfer function. The nB s subscript becomes mB R for the 
rotor-fir s t case 

(nBs)* circumferential Fourier mode of the impinging wake 
parameter used to describe the maximum deficit in the impinging wakes 
freestream or flight speed 
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U B> 

Uf(r*) 
- inc 

w v 

z 

z* 


^mid 

^lincr 


swirl "fireestream" between the stator and the rotor in a stator-first 
configuration 

swirl "fireestream" between the rotor and the stator in a rotor-first 
configuration 

relative ffeestream for the radiating blade section at i— 1* 

incident field's radial particle velocity over r=a + ; dimensional (Appendix A) 
axial coordinate for a cowl control point; nondimensional 

axial station of a typical radiating point for the incident field 

midpoint axial position for a liner ring patch (Appendix A) 

liner impedance (Appendix A) 


Greek 

a axially variable admittance of the liner, normalized by the fluid’s characteristic 

liner _ J 

admittance 1/pc 


Pax 2 I'M - 2 

Y m transform wavenumber for the radial direction, Eq. 5 c. 

Ap,AW unknown diffraction loading p'~p in the lifting problem, and unknown 

monopole strength w + - w' for the nonlifting problem; each of these two 
quantities takes on a v subscript when it becomes the circumferential mode of the 
parent variable; "+" and "-" denote r=a + and r=a 

7! liner loss factor (Appendix A) 

0 o off-axis directivity angle measured downward from the forward flight direction 

lifting kernel, or influence function linking the induced radial velocity to the 
virtual lift (Appendices A and B) 


iii 



.CAMBRIDGE ACOUSTICAL ASSOCIATES, INC. 


K/* 


p 

4 >. 4 >’ 

4»i( r+ ) 

Q 


kernel linking the radial upwash to the thickness source strength, or pressure to the 
virtual lift (Appendices A and B) 

kernel linking pressure to the cowl's dynamic thickness source strength 
(Appendices A and B) 

relative angle between the vortical wavefronts of the impinging wake's effective 
gusts and the chord surfaces of the cutting blades, following the nomenclature in 
Goldstein's book 12 


circumferential mode; short for mB R -nB s , the spinning mode index 


pitch angle for the B R rotor blade sections at radial station r* 
pitch angle for the B s stator blade sections at radial station r* 
background fluid density 

field and source values of the circumferential angle 
parameter to describe wake lean 
propeller's rotational speed 

frequency, ultimately mB R Q in the tonal problem investigated here 
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I ABSTRACT 

This project has had two specific technical objectives: (1) To develop a realistic three- 
dimensional model of tonal noise due to rotor/stator interaction, as the input field for predictions 
of diffraction and dissipation by a lined cowl; and (2) To determine whether the generator curve 
of that cowl, or duct, could be "steered" to yield substantially lower values of propulsor noise 
along the engine's fore and aft open sectors. The more general and important aim of the research 
is to provide the commercial aircraft industry with a useful predictive tool to help it meet its 
noise-reduction goals. 

The work has produced a tractable and yet realistic model of rotor/stator interaction noise. 
The blades in the fan stage are radially divergent, twisted, and of realistically wide chords to 
match the regime of high frequencies and speeds of the sound-production process. The resulting 
three-dimensional acoustic nearfield insonifies the interior wall of the diffracting cowl, whose 
shape, incidentally, does not affect fore or aft noise significantly. The predictions address both 
the fan's blade-passage frequency and its first overtone, "m=2". They indicate moderate acoustic 
shielding for the lined-cowl system; for example, see Fig. 17a. 
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H INTRODUCTION 

A. GENERAL 

The aims of the Phase I project have been twofold: (1) To produce, as an extension of the 
work in Ref. 1, a realistic three-dimensional theoretical model of rotor/stator noise from the fan 
stage of an Advanced Ducted Propulsor or Propeller (ADP 2 * 4 , Figs, la-c); and (2) To couple the 
resulting theory rigorously to a model that calculates the diffraction effect of a short cowl, or 
lined duct, on the propagation of that complex propulsor nearfield. 

This second goal was to be achieved through a coupling of formal unsteady lifting and 
nonlifting surfaces as mathematical replacements for the scattering unflanged cowl. Appendix 
A, excerpted from the Principal Investigator's Ref. 1, collects the integral equations that govern 
aeroacoustic diffraction for a straight duct of circular cross section. Myers and Lan 5 have 
recently implemented the time-domain version of the "lift" part of that theory, i.e., for an unlined 

cowl, in their own independent research effort. 

An important side question for the present project was whether the generator shape of the 
modeled cowl could be generalized theoretically, and used practically, to broaden the zone of 
null near-axis radiation of the lower order spinning modes resulting from rotor/stator 
interactions. Appendix B documents that analytical extension of the cowl's coupled lifting- and 
nonlifting-surface theories to include arbitrary warp in the duct s fore/aft direction. The cowl s 
cross section remains circular at every axial station. The speed/frequency regime is still Ref. 1 s 
high subsonic range for an arbitrarily noncompact cowl of "halfchord” L/2. 

The predictions of cowl-diffracted rotor/stator noise were to be provided for two kinds of 
observers whose perception of that sound would differ fundamentally for moderate to high 
values of the flight Mach number: (a) for passengers on the aircraft, relative to which the 
position of the noise-producing stators clearly does not change; and (b) for listeners on the 
ground, who would experience both the incident and the cowl-diffracted fields as sweeping by, 
i.e.. Doppler-shifted and with their rays amplified in the forward direction. 

The Phase I theoretical development has considered rotor/stator as well as stator/rotor 
interactions. The study just completed has resulted in four computer codes that predict cowl 
diffraction rigorously for noise from blades that are radially divergent, twisted, and with chords 
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diffraction rigorously for noise from blades that are radially divergent, twisted, and with chords 
that may be either aerodynamically compact or noncompact. Code 1 takes the stator to be 
upstream of the rotor, and the downstream radiating rotor blades are chordwise compact. Code 2 
is the same as 1 but for noncompact rotor chords. Codes 3 and 4 are the rotor-first, stator-second 
counterparts of codes 1 and 2 (codes 3 and 4 are of immediate interest to the civil aircraft 
industry, while codes 1 and 2 could have other applications). 

The side question regarding the potential effect of cowl camber or taper on on-axis engine 
noise has been answered, at least temporarily, through the development and exercise of offline 
canonical models rather than through the implementation of Appendix B’s exact theory for a 
doubly warped cowl in unsteady subsonic flow. That is because the canonical calculations 
concluded, relatively early on, that unrealistically high values of camber or taper would be 
required to achieve such off-axis kinematic "steering" of fan noise. The remaining Phase I 
research time was devoted to the thorough source modeling described in the previous paragraph, 
e.g., to the nontrivial analysis and numerical prediction of the nearfield velocities and pressures 
with which a set of radially divergent and twisted stator blades of wide chords would drive, 

aeroacoustically, a formally diffracting straight cowl. 

Readers not primarily concerned with theoretical developments should proceed directly to 

Sections V and VI. 

B. ADVANTAGES/FEATURES OF THE PRESENT LIFTING/NONLIFTING SURFACE 
APPROACH OVER ALTERNATE FLUID-FINITE ELEMENT MODELS 

(1) Our system of surface equations addresses one dimension less than the volume approach , 
with a corresponding potential gain in computational efficiency. 

(2) The lifting-surface theory developed under Phase I is of the acceleration potential type. It 
contains implicitly, and exactly, the shed and trailed vortex flows of the cowl in its role of ring 
wing. Those trailed and shed vortex systems do not spread for inviscid flow, i.e., they do not 
decay spatially with distance past the cowl’s outlet, or "trailing edge", so that their capture as part 
of the diffraction solution should be especially cumbersome via finite elements. Moreover, the 
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stator blades shed their own "incident" vortex system if the airloads resulting from the whipping 
effect of the rotor's wakes are properly modelled as the product of fully unsteady aerodynamics. 

Only by invoking q uas i-steady aerodynamics, which incidentally becomes untenable for 
the chordwise-noncompact regime of actual practice, could that pattern of directly shed vortices 
from the fan stage be neglected when calculating the volumetric flowfield. The present 
technique already accounts for these incident shed vortices, while the finite-element approach 
would again have to capture them through a laborious downstream gridding similar to that 
required past the cowl's trailing edge (the pole k+mBQL/2U m in the wavenumber 

transform of Ref. 1 's Eq. 57b contains analytically the field of induced vortex flows from the 
radiating blades). 

An important remaining problem in the overall research effort is the unified prediction of 
cowl diffraction and multiple scattering by all of the blades in the three-dimensional fan stage. 
That modelling extension would include the rotor’s forward transmission loss as a by-product. 
Blade scattering, or diffraction, is a fundamentally unsteady-aerodynamics phenomenon which, 
for the frequencies of interest and for the reasons just stated, would seem to be beyond the 
feasible computational capabilities of the currently used finite-element technique. That finite- 
element code would have to remain practically "non-unified" and therefore fundamentally 
iterative, or approximate (cf., for example, Meyer’s 7 nicely detailed paper on the importance of 
inlet reflections, calculated for one iteration). The volumetric approach would have to keep 
relying on an intermediate canonical code to furnish it with modal generalized forces that it 
would then "propagate" . The surface approach could eventually handle the problem of inlet and 
outlet reflections in a unified rather than iterative manner. That is because its influence functions 
may, and in fact already are, expressed spectrally both for the incident and diffracted fields. 

E a <-h such function appears in Appendix A and in Ref. 1 as a wavenumber transform in k . The 
combination exp (-ikz -zv4>) could serve in a formal superposition to determine the a priori- 
unknown incident nearfield at every diffracting blade strip in the fan stage. The resulting system 
of running "plane waves" incident on each such strip could be converted analytically into a 
diffraction load distribution for that airfoil section. That chordwise load distribution would 
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"return" to the three-dimensional medium as an elemental reflection and transmission whose 
strength would be found as part of the otherwise rigorous solution of the global 
cowl/centerbody/blades scattering problem. 

(3) The surface approach meets the radiation condition exactly through the theory s building- 

block Green's function; i.e., through the point-source solution of the converted wave equation for 
a uniform freestream U„ of Mach number The surface theory does not have to worry 

about spurious returns from infinity (although neither does the latest generation of volumetric 
techniques, which use infinite elements). 

(4) The lifting/nonlifting-surface solution of the cowl diffraction problem contains all of the 
inf ormation of duct acoustics, even though it is not formulated with a waveguide in mind. A 
post-processing of the combined incident plus diffracted fields inside the cowl would display all 
of the modal content of a leaky waveguide, with cut-on frequencies for each of the incident 
field's generalized forces, etc.; cf., for example. Fig. 4 from Ref. 8. 

(5) A more "buried" advantage of the surface approach is that its kernel functions have been 
analytically regularized, and that the singular parts have been found to be independent of the 
spinning mode number v = mB R - nB s {B R and B s are the number of rotor and stator blades 
in the fan stage). The submatrices associated with those singular subkemels may therefore be 
recycled from each v singular-equation system to the next. These singular subkemels show up as 
the first term on the right hand sides of Eqs. A.5 and A.7 . 

C. DISADVANTAGES OF THE PRESENT LIFTING/NONLIFTING-SURFACE 

APPROACH RELATIVE TO FINITE ELEMENTS 

Our approach is fundamentally boundary-element in nature. 

The main drawback of boundary-element numerical solutions to problems in unbound 
media is that the governing differential equation on which those techniques are based may not 
have arbitrarily variable coefficients. All boundary-element techniques rely on the analytical 
existence of a free-field Green’s function to convert the original volumetric differential statement 
into an integral one, expressed only over the domain’s bounding surfaces (to be sure, a Green s 
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function can usually be found approximately for the general case, e.g., by WKB 9 ). 

For the problem at hand, i.e., for the rigorous prediction of diffraction or scattering by the 
cowl, and eventually also for the coupled centerbody, the restriction of having a convected wave 
equation with constant coefficients causes the freestream and associated Mach number to be 
constant throughout the engine. The main immediate consequences of that are, or would be. 

(a) A centerbody that would not satisfy steady, flow tangency, i.e., that would remain 
hydrodynamically transparent to the axial steady freestream striking it The surface-only exact 
approach would lack the effects of ray bending which the missing steady streamlines would 
have contributed. 

(b) Swirl between the rotor and the stator generates a local convected wave equation with 
effectively variable coefficients. It is generally not possible to deflate it as usual to a non- 
convected wave equation with even derivatives through the usual changes in independent and 
dependent variables (or through the equivalent use of simultaneous Lorentz and Galilean 
transformations). The present approach thus cannot account for the propagation effects of swirl, 

at least not exactly. 

(c) The effect on acoustic propagation of a fan exhaust that is larger than the external flight 
speed U m is also not be easily handled by the present approach, although it may be possible, 
with some work, to adapt the surface approach to treat abrupt changes m JJ m from one constant 
value to another through the engine. 
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HI INCIDENT FIELD ANALYSIS 

A. CHORD WISE-COMP ACT AND CHORDWISE-NONCOMPACT SOURCE 

REPRESENTATIONS FOR TWISTED STATOR AND ROTOR BLADES 

A.1 General 

A major objective of the current study is to model mathematically not only the diffraction 
and dissipation effects of a short lined duct in a compressible freestream, but to generate also the 
incident field "p inc " driving that cowl from within. The analysis of this part of the problem will 
confine itself here to well-known mechanisms of tonal noise from turbomachines: (1) the 
chopping of stator wakes by a rotor, as displayed on the left half of Fig. 2a, and the whipping 
effect of that rotor's "viscous" wake on a downstream stator 10 , shown on the right-hand side of 
Fig. 2a n . This second type of unsteady-aerodynamics "interaction" is normally the only one of 
interest in the fan stage of current designs of civil-aircraft engines. Figs, la-c show modem 
examples without inlet guide vanes. 

The theory developed here for the cowl's-felt incident pressure covers both of the 

sources of noise in Fig. 2a, though one at a time as described in split form by Figs. 2b, c 12 . Our 
aim in considering both cases is to further commercialize the work by having it apply also to 
turbomachines not necessarily within the civil aircraft industry. The equations given below in 
Sections A.2 through B.3 address the less relevant case in Fig. 2b. A final section, B.4, provides 
the key for the transfer of symbols and parameter values for the conversion to Fig. 2c. 

A.2 Description of a Three-Dimensio n al Stator/Rotor Geometry for the Generation of the 

fowl's Internal Incident Field ('Stato r-First Rotor-Second Configuration} 

Reference 1 considered a single generic inflow inhomogeneity stemming from an axially 
aligned stator blade. The downstream rotor blades were multiple in number but lay flat on their 
plane of rotation; i.e., they were untwisted and perpendicular to the implied single upstream 
vane. Their sound-producing unsteady loads were the result of assuming a chordwise compact 
aeroacoustic field from root to tip even when actual turning and flight speeds, and chord sizes, 
mi ght have demanded otherwise. 

The geometry of the stator/rotor stage for the rotor-airloads model will now instead be as 
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depicted in cross section in Fig. 2b. There are B s stator blades rather than one. Moreover, they 
have an assigned pitch distribution u s (r * ) from root to tip, i.e., from r = R h to r = R t , to 
be able to treat later the rotor-first case of Figs. 2c and la-c (Section B.4 below). The rotor 
blades are now similarly twisted according to ), which, unlike ), will not be 

prescribable but rather will be deducible from the criterion that all blade sections in the rotor 
operate at zero angle of attack relative to the local freestream U r (r *)• All leading edges in the 
rotor will be unswept, viz., straight from root to tip relative to that local grazing freestream. 

The present study implements Fig. 2b to generate two separate models for the 

areoacoustic sources insonifying the cowl from within: 

(1) The kinematically elaborate extension of Ref. l's chordwise-compact rotor airloads 
to include the new geometric complexities for both sets of blades in Fig. 2b, 

(2) A new chordwise noncompact rotor-loads model that will again go 
geometrically with Fig. 2b. The rotor blades will be a set of "leading-edge" semi- 
infinite flat plates, each with a v R (r*) twist distribution from root to tip. The 
chordwise pressure distribution for each plate, at an arbitrary radial station r , will 
come from applying the "strip" assumption to that blade section. 

A main feature of Fig. 2b is to account for a new "freestream U^" just downstream of 
the stator. The unsteady-aerodynamics model for the source strengths will allow that meanflow 
to differ not only in direction with respect to the axial flight speed U^., as indicated in Fig. 2b, 
but also in terms of magnitude: U s> — U s> (r ) could be radially variable. 

This ends the qualitative description of the new airloads models of the rotor in its role of 
aeroacoustic source. The second half of the insonifying part of the theory is to "propagate the 
sound from these blade-based sources to the cowl's two-sided wall to drive the diffraction 
problem posed there. Fig. 2d shows the freestream model adopted during the current phase of 
work for this second part of the incident-field calculation: the swirl implied by U s> (r ) has 
been taken away and the uniform axial flight speed £/„ applies throughout The omission is 
typical of most "classical" solutions of the rotor/stator interaction problem; cf., for example, 
Envia's 10 Eq. 2.1 Id, where the substantial derivative for the convected wave equation contains 
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only the axial freestream component. A modem exception, though essentially two-dimensional 
or cascade in nature, is Hanson's recent analysis 4 of sound refraction through the tilted flow 
between the two sets of blades in our Fig. 2c. 

A.3 Mathematical Construction of the Impinging Wakes 

Fig. 3 displays the modeled spatial perturbation of the mean flow U s> (r) downstream 
of the stator blades. That "momentum-deficit" field will be Gaussian in shape and given by 


u(r,ty,Az) = u(r,*z)e 




( 1 ) 


where u (r, az) is the nonuniformity's maximum value and A (r, az) a free dimensionless 
function that regulates the circumferential sharpness of each modelled stator wake. Both u and 
A are typically taken to have a strong dependence on the axial distance az downstream of the 
stator disk. The value of az relevant to the present theory is obviously the axial gap between 
the rotor and the stator. Parameter (j) t (V) traces out the circumferential position of maximum 
deficit as a function of radius; i.e., 4> f (r) is the wake's "lean" angle as described, for example, 
by Envia. 10 

The Fourier-series expansion of u ( r , 4>, az) in the angular segment 0 < <f) < 2n/B s 
between adjacent stator blades yields a set of harmonic components u- given by 


u(r,<f>,Az) = Y, U n( r ’ AZ ) ein * 


(2a) 




_ 2 n/B s 

B 5 


w-(r,Az) = — f d$ «(r,$,Az) 

2u J 


(2b) 


A similar expansion again in (j> , but now over the whole circumferential range 
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0 < 4> < 2n , introduces the components u n (r, a z) for the counter "n" that are ultimately of 
interest in the airloads model (n without the overbar). For a stator with B g identical blades, that 
circumferential "n" spectrum filters down to a sparser nB ^ spectrum whose T , a z) 
components are related to the u- (r, az) values of Eqs. 2a,b by 


u n 



«jf (r,!.z)e 


-ire (nB s -n )/B s 


sin[7t(wB 5 - n )/i? 5 ] 
n(nB s -n)/B s 


( 3 ) 


The exponential and sine factors on the right side of Eq. 3 combine to give a Kronecker delta 
6 - when B-= 1 , as they should. 

TUI u 

The components normalized by the flight speed U m become 


= u (r) * e .„w s .* )/Ss sin[n(«B s -w)/BJ 

U m U„ 2 7 i jfe. 6 n(nB s -n)/B s 



-in$ -A(r)[fy - ^(r )] 2 

6^ 


( 4 ) 


The dependence of u and A on a z will be implied throughout. 

A.4 The Near Field Radiated bv the Rotor in Fig. 2a f or Both Chordwise-Compact 
and Noncompact Loadings 

We postpone stating the precise form of the rotor's source-strength solutions, and present 
now the generalization by the present theory of Ref. l's Eq. 57a: 





-imBfClt 

e * 



-i(mB R -nB s ) («|»-«|>*) 


PmB R -nB s ( a > Z ) ’ 


(5a) 


where the tilde on p inc on the right side, as in Ref. 1, denotes the harmonic spectral value of the 
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tonal signal, and the mB R - nB s subscript has turned out to be the circumferential modal 
counter relevant to the problem. This double spectrum p m B*- ti b s is* f° r i- e -» along the cowl s 
wall. 


^ B, , - m ) , . a ne ( r * *. £ ) s ™ [t*( r * )] 

= ' 7T / / dr — 


p uLl 


exp 


■ 

mB„QLM ' -- 

* . 

-I 

R m + k 

(z-z ) 


2c B 2 

J 



(5b) 


mB R -nB s y ' * m 


(y m O 


1 — mB p QLM ^ 
k + 


L \ 


2c P 


ax 


(mB R -nB s )cos[v R (r*)] 

sm [v R (r )] + ; 


where C is the contour shown in Fig. A. 1, where the product of the cowl's reduced frequency k 
and the flight Mach number become mB R ClL/2c. The new variables in Eq. 5b are 


/ mB D QL ^ 2 


Y m - 


N 


2c 


P 


pL* 2 


(5c) 


ax 


sin[p(r*)] = 


QLr*cos[u 5 (r*)] 

2 uy] 


(5d) 


sinlu^Cr’)] = 


- XJ S> (r * ) sin [u 5 (r * )] 

Zr 


(5e) 


uy) 
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cos[u^(r*)] 


^>( r *) C0S K( r *)3 

uy) 


(50 


uy > 


N 


t/i(r*)cos 2 [u Jr*)] + 



U s y)sm[v s (r*)] 


2 


(5g) 


Variable r* is nondimensional throughout Eqs. 5b-g; its original dimensional form has been 

normalized by the cowl's axial half length L/2. 

Eq. 5 b returns to Ref 1 'sEq. 57afor v R (r*) = *-12, u 5 (r*) = 0,and U s> = 0 in 
Eq. 5 d (so that U r (r *) = Qr *L!2 and sinful (r )] = 1 ). The second term within curly 
brackets in the it integrand accounts for the new component in the circumferential direction 4 > 
for the unsteady blade loads. It plays the role of the "drag" dipoles in Ref. 12's Eq. 3.1 19. 

The group {u^r * >7^ (r *, k ) sin [p(r * )] }/p U* comes from the aerodynamic 
theories of the chordwise-compact and noncompact models. The k variable will drop out of the 
argument list of the aerodynamic transfer function in the former case. 

A 5 Chordwise-compact and Noncompact Airlo a d Solutions for the Rotor Blade? 

A.5.a The Improved Compact Model 

Dittmar et al. a have investigated the potentially beneficial effects of stator 
systems that were designed to be drastically chordwise noncompact. But that study will turn out 
to be qualitatively relevant also to more conventional modem stator designs, essentially because 
the aeroacoustic parameter c»)c/ 2 c of these recent stators sometimes approaches the 

effectively infinite value of those in Ref. 2b. 

We report here only final results. Beginning with the chordwise-compact case, the 

normalized unsteady sectional lift on the rotor blades for R h < r <R f is 
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P U±L 


— i • 2^ • r* cos[d (r*)] 

4 J t R t 


• S' 


^ «5 5 QZ,c(r*)^ 
4tf r (r*) 



- n)/B s sm[n{nB s -n)/B s ] 
n(nB s -n)/B s 
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-ink-Afr-Wtf 
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where c(r * ) is the blade chord at the radial station r * normalized by L/2, and S is the Sears 
function. J ( is short for the advance ratio n UJiQLRJ 2), with R { similarly normalized by 
the cowl’s halflength L/2. 

A.5.b The New Analytical Non compact Model 

Variable z * in the exponential term in Eq. 5b marks the axial positions of the 
source distribution. For the compact case z has a single value so long as the propeller s set of 
quarter-chord lines has no rake, i.e., so long as it is coplanar. That was implied in Ref. 1 s Eq. 
57a. Variable 4>* in Eq. 5a above plays a similar role for the circumferential direction: 4>* = 0 
refers to a circ um ferential datum position for the points on the surface of each rotor blade. <|> is 
dispensable when the loading on the rotor blades is chordwise compact and when the blades are 
straight rather than swept; i.e., t|)* is to sweep what z * is to rake, for a set of chordwise 
compact rotor blades. Ref. 1 considered unswept blades only and therefore set (j) to zero 
throughout its analysis. 

Both z * and <|)* become necessary in Eqs. 5a and b when the blade airloads are 
chordwise noncompact even if the blades are straight and their rake is zero. One begins by 
setting z * = z* t + z t * 0 , where z* t is the "global" axial position of a chosen point in the radiating 
rotor, and z* i0 measures the local axial positions of blade points away from there. Variable (j) is 
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now understood to be similarly local without need of the "Co" designating subscript Then the 
relationship between z *, <J>* and the chordwise variable X * , which appears marked on the 
uppermost rotor blade in Fig. 2a, is 


-r*4>* = x*sin[i>*(0] 


(7a) 


Zgi = z lo = ^*cos[u^(r*)] 

And the phase factors in Eqs. 5a and b become 


(7b) 


exp 


-i(mB R -nB s )( 4>-4> ) - 


imB r QL M m (z-z * ) - 


ik(z - z * ) 


= exp 


-i(mB R -nB s )i f> 


2cP« 

imB R Q L (z - z g[ ) - 


2c P 


- ik(z-z gl ) 


(8a) 

/%(**) , 


ax 


where 


Ph(x*) = exp ^ x* 


(mB R -nB s ) 


sm[i>»(r*)] 


- mB p QLM 
k + 2 S. 

2cpL 


cosfu^Cr*)] 


(8b) 


The group w - T- sinp/ p U^L in Eq. 5b's r * integrand now changes to 


% nB s *nB s 
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J 
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Ph(x *) 


( 9 ) 


Q uant ity Ap(x* denotes the loading on the blade section in question normalized by the 

effective gust amplitude . This ratio is given analytically by Eq. 37 of Ref. 13 or by Eq. 15 

£ 

of Ref. 14. In the former case the sweep angle is set to zero with the interpretation of the 
formulas in both studies considering a gust parallel to the leading edge. 

We ghall further simplify those canonical results by stripping them of their trailing-edge 
contribution, given that dial contribution acts as a small correction to the dominating leading- 
edge part for the high-noncompact regime of interest here. Based on the 
exp (- (= exp (- imB R Q t)) temporal behavior adopted in Eq. 5a, then 


A p£) = 


iic/4+i- 


nB s QLx 

2c(l+A^) 
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- 2p U r 


nB c 



nB~Q.L . 


5 X 

2U r 


( 10 ) 


Substitution into Eq. 9 yields, with Eq. 5d, 


U nB'(r ‘ ) T nB.( r ’> *) shl D* ( r ’ >1 U nB^ r * > Q L 
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where 


m = 


nB s QL 

2c(l+M r ) 


- mB p QLM 
k + — ^ ^ 


2cp 


or 


cost^Cr*)] 


) 


( mB R -nB s ) 

sinfu^r )] . 


E is a Fresnel integral defined by 


(12) 
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£(<10 = / 
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(13a) 


For/(Jfc) real and less than zero, the following would take over in Eq. 1 1: 


E{c{r')f(k)} E'{-c(r')f(k)) n3w 

; — ' j 

vTw 'i -/(*) 

where the asterisk on E * stands for complex conjugate. The final step for the present 
development is to bring Eq. 4's expression for !U m into the first term on the right side of 
Eq. 11. We omit it for the sake of brevity. 

B. FARFIELD EXPRESSIONS FOR ^ 

B.l General 

We shall follow Ref. 1 by considering again not one, but two separate farfield 
calculations that produce clearly different results when the flight Mach number is high 
subsonic (otherwise the two farfields coincide): (1) the directivity pattern perceived by observers 
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traveling with the aircraft, and (2) the corresponding field for ground-based observers. Farfield 

(1) is the one of interest for studies of transmission of engine noise into the cabin while in flight 
(assuming the distance from the engine to the cabin is large in terms of wavelengths). Farfield 

(2) becomes relevant for predictions of the environmental impact of the ducted or unducted 
propulsor on listeners on the ground, e.g., dwellers near an airport. 

The following two subsections document the final results of both types of field for the 
chordwise-compact and noncompact models of loading on the radiating rotor blades. Those final 
results for ^ are interesting on their own right, but our main objective for generating 
them here is to complete the field of total pressures, incident plus diffracted, which are expected 
to have low levels in the shadow of the cowl's wall. 

B.2 Farfield for Incident Rotor Sources in th e Aircraft's Frame of , Reference 

Ref. l’s Eqs. 61a, b defined the directivity ang le 0 O as tan' ! (- P^Wz) and the 
nondimensional range R o to the observer as sfr^+z 2 / P^ ■ Neither variable was therefore 
completely physical in that both invoked P^ for convenience. We redefine them here to be the 
more directly useful quantities tan '\-r/z) and \/r 2 +z 2 , respectively. E.g., R 0 is now given 
by Ref. l’s symbol for R o divided by p^l-M^si n 2 0 o , a new factor taken now into the 
denominator of the right-hand sides of Ref. l's expressions for incident and diffracted fields. 

These name changes also have the effect of turning Ref. l's " COsQ^ " into 

cos 0 O / ^1-M^sin 2 0 o , and Ref. l’s ’’ sin0 o " into P^sinO,, /^l-M^sin 2 ©^ The final 
result for now is 
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mB R QL 

4 

2c 



[“.('•')] + " Bs> cos [«*(»• *)] 


The aerodynamic transfer function T ^ ( r , 0 O ) for the compact case ceases to be a function of 

s • 2 . 

0 , and Eq. 6, as is, supplies then the expression for ^hB 5 P • ^ or 

noncompact case, one finds that this group function now takes on the form 
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B.3 Incident Farfields in the Ground Reference Fra me 
Instead of Eq. 14, one now finds that 
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Eq. 6 a gain provides without changes the r -integrand factor '^hB s ^ ^ or ro * or 
blades with chordwise-compact loads. For the noncompact case, one obtains 
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where now 


AV = 


nB~QL 
± + 

2c(l+M r ) 


mB R QL 

2c 


costu^Cr*)] cos0 o 


(mB R -nB s ) 
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sin [ 1 ^( 7 * *)] 


(17b) 


B.4 Gon version of the Blade-Source Analysis in Sections A.7. - B-3 to Apply tp the 
Rotor-First Stage in Fig. 2c 

Again, the analysis of the fan stage as a source of noise, described above in Sections 
A.2 - B.3, addressed Fig. 2b's case of inlet guide vanes preceding a rotor. The purpose of this 
section is to modify all of that work to make it also fit the rotor-first arrangements shown m cut- 
away form in Fig. 2c and for whole engines in Figs, la - c. This second configuration has been 
the subject of numerous recent analyses, cf. Envia 10 for a "classical” problem of the Wiener-Hopf 
type, and Kousen & Verdon 15 for a CFD cascade solution in two dimensions. 

In terms of parameter values, the rotor-first arrangement differs from the stator-first as 

follows: 

(a) The rotor's viscous wakes clearly belong to a higher Reynolds number flow, and their 
"laminar" widths are accordingly narrower; 

(b) Their "lean" could be greater than for the stator-first case (cf. Envia 10 ), because ... 

(c) Both the axial and circumferential components of the rotor's exhaust freestream 
U^if) would be numerically greater than U s> (r ) for the stator-first case: Momentum 
conservation requires that the rotor's performance thrust and torque be reflected in 
"kicked-up" values for these quantities; 

(d) The final exhaust angle for the stage would now have a practically small chordwise- 
average value, unlike the values suggested in Figs. 2a and 2c, and more like those in 
Fig. lc. 

In terms of the character of the unsteady-aerodynamics problem, and of the resulting 
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noise from the downstream blades, the stator-first and rotor-first pictures differ as follows: 

(a) In Section A.3 the stator's wakes were spatially stationary with respect to the 
downstream rotor, which cut through them. Now they are obviously in the nature of convected 
gusts passing over the downstream stator blades at the collectively sampled blade-passage 
frequency plus harmonics, which remain the same as on Sections A.2 - B.3; i.e., given by 
jnB D Q - The airfoil-felt disturbance wavenumbers are now given by mB R Q/U r rather than 

by nB s Cl !U r . 

(b) The character of the convected wave equation in the "propagation" part of the 
problem does not change, because the formulation in Eq. 5b neglects the swirl component of 
U and U D for the two situations in Figs. 2b, c, and takes as applying throughout. I.e., 
as modelled, for example, by Envia 10 (again, his Eq. 2.1 Id), and in Fig. 2d here, the blade- 
radiation part of the analysis is the same for rotor-first and stator-first configurations. 

One concludes from all of the above that the following conversion key carries the 
stator-first source theory developed in A.3 - B.3 to apply to the rotor-first configuration. 

(1) Exchange the symbol B R , which appears in the numerator of the coefficient 
immediately to the right of Eqs. 14 and 16, for JB^the number of now-radiating stator 

blades. Also switch in those expressions the local pitch angle ^(r) of the radiating 

blade section to be u s (r ) ; 

(2) The index "nB s " of the group that denotes the submodel for the unsteady 

load on the radiating blades now becomes mB R , the index for the new convected 
wavenumbers striking the downstream stator; 

(3) Accordingly, in the compact-loads model of Eq. 6, the nB s part of the Sears^ 
function, and B s itself, become mB R and B R , respectively. The blade chords c(r') 
now refer to the stator rather than to the rotor, and again v s (r )~+ \) R (r ) ; 

(4) For the noncompact load results in Eqs. 15a,b and 17a,b, items (2),(3) just above 

apply as well in Eqs. 15b, 17d. 

One final remark: Sections A.l through B.3 dealt with a radiating rotor with blades that 
obviously did not reach the cowl’s inner wall. I.e., the tip radius R , was less than the cowl's 
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radius "a". The same is clearly not true for radiating stator blades, which <k> connect structurally 
to the cowl. For reasons of numerical convergence, however, we assume here that the effective 
radial extent of stator radiation is again less than "a". The justification is a consistent use of strip 
theory for the now-downstream stator blades, together with a picture of rotor wakes that do not 
extend radially outboard to the cowl's inner wall at r—a. 

IV A CANONICAL MODEL OF COMBINED LEADING- AND TRAILING-EDGE 

DIFFRACTION FOR THE COWL: ACOUSTIC SHIELDING BY AN "UNWRAPPED" 

DUCT 

A. MOTIVATION 

Figs. 4a,b show the geometry under consideration: Fig. la's cylindrical cowl has been 
unwrapped into a strip of width equal to the original cowl chord L. The length of that strip, 
which corresponds to the cowl's original spanwise extent 27taas a ring wing, will now become 
irrelevant as per the rules of the Geometric Theory of Diffraction (GTD) 16 , which this section 
adopts. GTD concludes that only one point on the cowl's inlet and outlet circular edges plays a 
role in the diffraction process: That which is closest to both the interior insonifying source and 
file external receiving field position. 

This somewhat narrow application, or interpretation of that theory , obviously neglects the 
multiple reflections that the cowl's actual waveguide geometry generates in its interior (more on 
this under Section V.B. Id). The unwrapped shape in Fig.4b provides a spurious infinite-medium 
reservoir "below" the source. The true geometry in Fig. 4a instead confines the acoustic energy 
of those lost "lower" rays to the cowl's finite internal cylindrical space. The duct s interior 
becomes a classical "whispering gallery" to those rays, which eventually work their way to the 
two "p rimar y" edges shown in Fig. 4b and affect the indicated listener through additional edge 

diffractions and possible specular reflections. 

The unwrapped-cowl model is nonetheless useful because it provides an estimate of the 
best shielding performance possible under the assumption that the duct's waveguide effects, i.e., 
its multiple reflections, etc., could be completely eliminated through absorption by a liner. The 
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r emainin g physical picture would ideally then justify the unwrapped-geometry limit in Fig. 4b. 
These long overdue canonical calculations will also help interpret the formal predictions of duct 
diffraction in Section V. 

The insonifying field in this canonical sub-study will be that of a single three-dimensional 
monopole source of volume- velocity strength Q 0 . The same fxeestream of subsonic Mach 
number of Figs. 2b-d will continue to tilt all of the acoustic fields in the downstream 
direction: incident, reflected, and edge-diffracted (more on ray tilting in V.A). 

B. BRIEF BACKGROUND OF SCREEN-DIFFRACTION PROBLEMS WITH AND 

WITHOUT FLOW 

The classical problem of diffraction by a rigid screen in a still medium produces a field 
that applies eq uall y to the velocity potential and the pressure. That is because the two quantities 
differ then only by a multiplicative constant, namely, I to p , as in the statement 
fap = _ (A O) /dt = ZO) p A 0 , which relates the "diffraction load" A p across the screen 
to a similar jump in potential <£ . Either dependent variable is furthermore found to vanish as 
the square root of the distance to the screen’s edges. And both p and 3> are zero over the 
coplanar extension of the screen. 

The presence of a mean flow changes that picture radically. A leading-edge screen then 
generates a diffraction load A p that blows up on approach of the screen's vertex, while A $ 
acts just as it did in the flowless problem (Fig. 4c). The reason for A p 's change in behavior is 
that now the relationship connecting it to A $ contains a spatial derivative: 

A/? = - pD(A<&)/Dt, which on approach of the edge at z=0 yields 
A p ~ - pL^dCASO/ dz ~ z A/2 , since A$ - z 1/2 as before. 

If the screen is a "trailing-edge" one, with M^. flowing past it as shown in Fig. 4d, then 
A p is zero along the screen's wake extending to Z = +°°, again as in the flowless problem 
(where there was no wake), while A $ should now be finite rather than zero there. 

Candel 17 , and then Amiet 18 , exploited these well-known features, of flow versus nonflow 
problems for the screen geometry, to construct leading- and trailing-edge pressure solutions for 
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M * 0 for incident fields of the "appropriate" kind. The more relevant of the two studies is 
ax 

A mi Kt's because it addresses the point-source problem in our Fig. 4a (Candel considers a plane 
wave incident from infinity). 

By incident fields of the "appropriate" kind we mean the following: The driving point 
source in Amiet's leading-edge solution could afford to be a standard velocity source of volume 
strength Q 0 . The p rimar y dependent variable in the boundary value problem was accordingly 
the velocity potential, since its particular solution (without the barrier) could be immediately 
connected to the point velocity source through a gradient operator that did not include the 

freestream C/^. Amiet then obtained, or could have obtained, the associated pressure field 

through a standard post-application of the operator - p D/ Dt to the independent time and axial 
field variables. Summarizing then, the leading-edge problem for 0 has a straight-forward 
solution for the velocity potential and for the pressure field, for a standard velocity source of 

M 

strength Q q generating the incident field. 

The same is not true of the trailing-edge screen in Fig. 4d here. There is then no analogy 
with a flowless problem unless the primary dependent variable is the pressure throughout the 
picture. That includes an exchange of the true velocity source driving the system, Q 0 , for an 
artificial "pressure" source, say of strength " P 0 ". Amiet's (and Candel s) leading- and trailing- 
edge solutions are therefore irreconcilable because their source fields are incompatible. 

C. POSTULATED RECONCILIATION OF EXISTING LEADING- AND TRAILING- 
EDGE CANONICAL SOLUTIONS 

One replaces the pressure source P Q cited above with the product of - p Q 0 and the 
"source-adjoint" version of the standard substantial derivative: 


PrT-PQc 


'±-1 J ^ 

v dt Sz* 


(18) 


25 



CAMBRIDGE ACOUSTICAL ASSOCIATES, INC 


The spatial derivative is with respect to the source position z* , and the coefficient in the 

ordinary substantial derivative has been replaced by - U m . 

The application of the right side of Eq. 1 8 to Amiet s trailing-edge solution for ^ 0 , 
after removal of his implied symbol 1 P g , contains the following features and we therefore argue 
is the correct solution for the trailing-edge problem driven by a velocity source: 

(1) The problem's particular solution is recoverable as the field-point vector X approaches the 
incident source point x' , given that d/dz * = -d/dz . Le, the solution captures radiation by 
the velocity source in free space, a weak but necessary requirement; 

(2) The resulting solution still satisfies the Kutta condition and flow tangency. The new 
construction has simply exercised its right to interpret the original complete solution as the 
barrier problem’s fundamental influence function. The operation in Eq. 18 builds up a more 
complicated field from that original total "Green’s function"; 

(3) A $ con tin ues to be nonzero along the screen's wake. The inversion of the standard 
substantial derivative -pD/Dt, through the integral operation 


ioz/U^ 2 

- — 1 f 


dz 


-toil U m 

p uT 


9 


does not "undo" the right side of Eq. 18 when z refers to the axial positions of the virtual 
sources that make up the screen's surface. Eq. 18 does not "trivialize" the associated solution for 

the potential: the shed wake stays (recall Eqs. 13-17 of Ref. 1). 

We avoid writing down the complicated result of applying Eq. 18 to Amiet s trailing-edge 
solution 1 *. The differential operation yields derivatives of Fresnel integrals, etc., that provide no 
great new insights into the nature of the diffraction process. Section E below will present sample 
results from the numerical implementation of that unstated expression. 
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D. HIGH-FREQUENCY JUSTIFICATION OF A COMBINED LEADING- AND 
TRAJLING-EDGE DIFFRACTION MODEL 

The end products of Section C are leading- and trailing-edge solutions for * 0 
driven by the same velocity source of strength Q 0 • They may now be referred to a common 
coordinate system, eg., to one centered at the field position of the driving monopole source, and 
added to generate an aeroacoustic diffraction model for the wing of finite chord earlier discussed 
in Fig. 4a. The implied assumption is that the edges decouple aerodynamically for the relatively 
high frequencies of interest, viz., for co L! (2 c (3^) in the range of 20. Cf. Peake, 19 for a 
recent development of a leading-edge solution for a quarter-wing in unsteady supersonic flow, 
based on Landahl's 20 chordwise decoupling of edge behaviors. Ring-shaped signals created at 
our cowl’s leading edge directly by Q 0 would also spread downstream quasi two-dimensionally, 
and their amplitude at the trailing edge would have decayed at least by the factor 
{coZ/[c(l+M )] }^ 2 . Pressure signals similarly generated at the trailing edge would be 
cut down by the even smaller multiplier {( oL/ [c(l -M^)] }' ,/2 on reaching the leading edge. 

E. NUMERICAL PREDICTIONS OF NOISE SHIELDING BY THE , TWO-EDGED 
CANONICAL MODEL FOR * 0 : EFFECTS OF COWL CAMBER AND 
TAPER 

The set of canonical diffraction calculations about to be discussed generalizes those of 
Dittmar 2 *, who also invoked diffraction by a rigid screen in the context of the cowl but 
considered only the flowless case. 

Fig. 5a shows the directivity pattern of incident (p inc ) , and inddent-plus-diflfracted 
(total) pressures for a leading-edge screen. This solution is immediately available from Amiet's 
work. 1 * The three- dimensi onal point monopole is 30° off the screen at a dimensional distance 
R* away from the vertex such that (x)R lc 10. The picture is the same as the sketch on the 
right side of Fig. 4b but with the trailing-edge point set at z = +°° rather than at L/2. is . 8 . 
The upper part of Fig. 5a's directivity pattern contains the shadow zone. The lower part contains 
the spatial region of specular reflections and a transitional zone connecting it to the shadow. The 
directivity pattern of the incident field would be a perfect circle were it not for the freestream, 
which distorts it into an oval shape. The gap between p mc and the total pressure on the shadow- 
zone (upper) half of the figure indicates the degree of shielding provided by the rigid barrier as a 
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zone (upper) half of the figure indicates the degree of shielding provided by the rigid barrier as a 
rigid leading-edge wing. 

Fig. 5b plots similar directivity patterns for the trailing-edge solution, with the same 
velocity source Q 0 now supplying the incident field for the first time (cf. the analytic 
construction under Section C above). 

Fig. 6 adds coherently the two total-pressure fields in Figs. 5a,b and thereby arrives at the 
diffracted field of the wing of finite chord L displayed in Fig. 4b. Whereas the shadow-zone 
patterns of the individual leading- and trailing-edge cases in Figs. 5a,b were smooth, their 
combination in Fig. 6 now contains the fringe typical of two effective "diffraction sources", one 

at each of the wing's two edges a distance L apart. 

Fig. 7a checks that Fig. 6 contains the right number of lobes within a sector roughly 
defined about the normal direction in the shadow zone. Fig. 7a plots the combined phase factors 
of two sources of arbitrary but equal amplitude, again separated by L. Fig. 7b checks the parts of 
Figs. 5a,b that are less physically relevant to our purposes though important as diagnostics: the 
roughly-equal lobe patterns of specular reflections generated by both semi-infinite screen 
problems on their illuminated side. One of the two curves in Fig. 7b is the total pressure from 
the leading-edge solution. The other is from the even more primitive problem of a source near an 
infinite wall, with = .8 again. The lobe patterns differ at the beginning of the zone of full 

illumination, as they should, but then quickly merge as required. 

Fig. 8 generalizes the two-barrier construction of Fig. 6 by staggering the two screens to 
model the diffracting ends of the generator curve of a "truer" geometry, shown below as a bell- 
shaped cowl of the same axial length as the straight duct in Fig. 6. The slope of its generator is 
zero at both the inlet and outlet. The source is now closer to the trailing edge; i.e., 

(i>R ’/c — 3.34 . The fringe pattern is less textured because the diffracting strengths of the two 
ends differ substantially from what they were in Fig. 6. The overall depth of the shadow, 
however, is about the same as in Fig. 6. 

Fig. 9 "pushes" the two-edge screen model by tilting the edges relative to each other. 
Otherwise the picture is the same as in Fig. 8. The tilted screens are a "tougher sell" now as a 
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model for the cowl shape, as depicted at the bottom of Fig. 9, because they require that the 
ffeestream also be tilted, but in different directions following the edges. The shadow softens, 
especially near the horizontal direction, because the source is now more "exposed”: the leading 
edge's angle relative to the source is now greater by 20° than it was in Fig. 8; the trailing-edge 
angle is greater by 25°. 

The main conclusion of Figs. 6 through 9 is that over a fairly large sector about the 
horizontal none of the three two-edge models provides substantial shielding of noise. The 
stipulated use of "steering" of the cowl's virtual dipoles to reduce on-axis noise is therefore not 
really practical unless the fundamental mode index B R - B s of the incident field is very high. 
The conclusion holds even for the relatively high realistic frequencies investigated through these 
canonical models. Fig. 10 summarizes why this is so. The remaining calculations will use the 
rigorous cowl-diffraction model to explore other means of achieving reductions of near-axis 
noise from ducted rotor/stator interactions, e.g., through liner dissipation coupled formally to the 

short cowl's diffraction effects. 
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V SAMPLE PREDICTIONS OF FORMAL COWL DIFFRACTION FOR STATOR- 
FIRST AND ROTOR-FIRST SYSTEMS 

A. KINEMATICS OF RAY TILTING IN A UNIFORM FREESTREAM 

The purpose of this subsection is to add a few remarks to the canonical groundwork of 
IV. The new co mm ents will help interpret the calculations discussed below for the rigorous 
diffraction of fan noise by the cowl. The issue at hand here is the expected ray-tilting effect of 
the uniform freestream on the cowl's pattern of reflection and diffraction lobes. 

Fig. 11a shows a pair of typical rays shooting out of an arbitrary point of radiation in the 
fan stage. One ray travels upstream and the other downstream. Both do so at the standard speed 
of sound "c" relative to the medium that convects them. The figure suffers no loss of generality 
in as signin g the cowl's midaxis station to the two rays' point of origin. The cowl s inlet ring 
appears as a point, at a "visual", or nominal angle 0 ™ nom with respect to the source. The 
corresponding outlet value is 0*%,; it is equal to the inlet value for the special arrangement 

considered in the figure. 

The lower part of the picture gives the solutions for the ray tilt needed to clear the cowl in 
both upstream and downstream directions. The upper part of the figure displays the cowl's 
optical, or geometric shadow Aef 1 ^ along with a sector of equal width spanning the tilted 
ray directions. To be precise: the cowl's acoustic shadow, in the presence of the flow, remains 
symmetric about the midaxis and matches the optical shadow. However, within the upstream 

J P . p inlet 

complement of the figure's uppermost sector, i.e., for field angles 0 O satisfying U Q 0% nom , 
the directional content of the unblocked ray field would be limited to the angle 0 

A corresponding observation applies to Fig. 1 la’s downstream direction, reconsidered in 
Fig. 1 lb for a typical ray that does nol clear the cowl. The contents of ray directions for incident 
and cowl-reflected rays include the indicated sector, which due to the flow is greater than the 

nominal. 
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B. RIGOROUS PREDICTIONS OF COWL DIFFRACTION FOR SECTION HLB'S 

INCIDENT FIELDS 

B.l Calculations for a rigid -walled (unlined) Oral 
B.la Parameter Values 

The parameters for the cowl and the fan stage are: The number of rotor blades B R is 16. 
The number of stator blades B s is 22 (following the Advanced Ducted Propeller case in Meyer’s 
paper 7 ). The fundamental value of the spin-mode index V = mB R - nB s is therefore -6, for the 
first temporal harmonic counter m=l and the stator's first Fourier counter n=l . The flight Mach 
number M^ will be .4 for all of the calculations. The rotor’s RPM is 4000, yielding a Blade- 
Passage Frequency (BPF) of 1067 Hz. The cowl's length L is 8 feet and its dimensional radius 
aL/2 is 2.5 feet. The resulting width A 6 f 4 " of the cowl's geometric shadow is 78° for a 
blade-tip radiating point (see engine insert at the top of Fig. 12). The calculations do not include 
the engine’s afterbody, which is transparent both hydrodynamically and acoustically. The value 
of(o>/cXaL/2) is 15.2. Wake lean is neglected for simplicity's sake: 4> f (r ) — 0. The wake 
parameters "A" and ul U m have the constant values (2/it) 2 log e 10 and .2, respectively. 

The present lifting- and nonlifting-surface theory of cowl diffraction does not invoke duct 
modal concepts explicitly, though it nonetheless "contains" them. I.e., the new theory makes no 
explicit use of the fact that at the BPF at least three radial modes have cut on in a hard-walled 
duct lacking a centerbody: the eigenvalue problem for an infinite rigid-walled cowl is 
j;( y t d) = 0 , which for v=*6 yields Yl a=7.2, 1 1 .6, 15.4, .... for l=l,2,3,... The associated cut- 
on values for (o/c)(aL/2) are given by P„Y,a (see, for example. Fleeter 3 ), which for M ax =.4 
(P»=.92) gives that (co/c)aL/2=6.6, 10.7, and 14.2 as the three relevant values below the 

fundamental operational value of (B R Q /c)aL/2, i.e., below 15.2. 

Regarding chordwise compactness versus noncompactness for the radiating blades: They 

will be chordwise compact at the BPF and noncompact at the first overtone, 2xBPF. This 
modelling decision is based on a conclusion reached by Amiet 21 for an airfoil's unsteady lift at 
high frequencies and speeds. It turns out that values of (to /c) [c (r) L/4]/(l-M r ) greater 
than n/4 justify the use of the leading-edge solution alone, which in the present case is the 
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chordwise pressure distribution given in Eq. 10. Values of the parameter lower than ic/4 
conversely call either for a multiply corrected noncompact solution, or for a fallback on the 
compact Sears solution given here in Eq. 6. The chord dimensions used in our calculations are 
such that the BPF value of Amiet's chord parameter is roughly .54 for airfoil sections near the tip 
of the radiating blades (i.e., lower than n /4 ). The appropriate model is therefore a chordwise- 
compact one. At the first overtone, however, the parameter value doubles to 1 .085 and crosses to 
the noncompactness regime of Section HI.A.5b. 

B.lb Predictions for an Academic Stator-First, Rotpr-SecondFan Stag e 

Fig. 12 documents an example calculation of diffraction by a rigid cowl for the study s 
only reported stator-first case. The frequency is blade-passage. The reference frame of the 
directivity pattern is that relevant to community noise (Section m.B.3). The solid curve displays 
the rotor's incident field in the absence of the cowl (201og, 0 of Eq. 16, etc.). The dashed curve is 
the total field obtained from the phased sum of incident and cowl-diffracted fields. The shielding 
performance of the rigid-walled cowl is mixed, with the original incident field intertwining the 

"total" curve. 

B.lc Predictions of Diffraction of Stator Noise by an U nlink Rigid Cpwl; m=l, the BE E 

The remaining calculations address the rotor-first, stator-second configuration of current 
interest to civil aviation. The flow parameters have been adjusted artificially to keep Fig. 12 s 
radial distribution of rotor-blade twist for the now-radiating stator blades. The airfoil disturbance 
wavenumbers are now given by mB R Q/U r (0 rather than by nB s Q/U r (r ), as explained in III.B.4. 
Again, U t (r*) is artificially the same as in Fig. 12, just as the shape and amplitude of the wake, 

etc., are kept the same. 

Fig. 13a plots incident and total fields for a ground observer. The frequency is again 
blade-passage and the radiating stator blades are chordwise compact. The cowl's mixed shielding 
performance is similar to that of Fig. 12. The lower overall levels of this rotor-first calculation 
relative to those in Fig. 12 are chiefly attributable to the difference in values of the relevant wake 
harmonic u n versus u R , which in Fig. 12 had a much higher value due to the fact that B R = 16 

is significantly lower than B s = 22. 
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Fig. 13b compares Fig. 13a's ground prediction (dashed) for the total field to that on the 
aircraft’s frame of reference (solid). Even for the relatively low Mach number of these cases, i.e., 
for M iX =.4, the two patterns differ considerably. Community noise is significantly higher in the 
flight direction while cabin noise gets washed downstream, both according to expectations. Fig. 

13c shows similar results for the incident field alone. 

Fig. 13d is a diagnostic comparison of the cowl's complete scattered field and the stator 
blades’ incident field. The figure has tentatively labelled the parts of the directivity pattern 
associated with the simple, one-bounce reflection phenomenon depicted in Fig. 1 lb and 
originally explored quantitatively in Figs. 5a,b and 7a,b. The number of lobes in the specular 
zones makes rough sense upon taking account of the difference in frequencies and Mach 
numbers between Fig. 13d and the earlier canonical results. The purpose of Fig. 13d is both to 
support and question the rigorous total field prediction in Fig. 13a: The scattered field appears to 
have the correct spatial structure, and yet the predicted shielding in Fig. 13a is far from the 
impressive conclusions reached by the canonical model. A hypothesis for the rigorously 
calculated cowl performance will be offered below for the unlined case. 

B.ld Predictions of Cowl Diffraction of Stator N oise bv an Unlined Cowl; m=2, the Fi rst 

Overtone of the BPF 

Figs. 14a-e are the m=2 versions of Figs. 13a-d: The frequency is now twice blade- 
passage and the radiating stator blades have become chordwise noncompact from root to tip, 
thereby calling for a computational switch from Eq. 6 to 1 1, etc. The spin-mode index v has 
climbed from -6 to +10 for V’ still equal to 1. Fig. 14a shows incident (solid) and total (broken) 
fields for an observer on the aircraft's reference frame. Fig. 14b is for a ground-based listener. 
The levels in both figures are considerably lower than those in Figs. 13a,b for the BPF, for this 
particular sample propulsor case. 

The cowl's shadowing performance displayed by Fig. 14a is somewhat better than that of 
blade-passage. It is still, however, below the expectations raised by the canonical model. Figs. 
14c,d compare again the incident to the total field for aircraft and ground-based observers. Fig. 
14e contains the overtone's diagnostic plot of incident and cowl-scattered fields. The upstream 
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and downstream lobe patterns for specular reflections are now clearer, and therefore more 
credible than they were at blade-passage. The potentially "shadow-producing" central zone is 
likewise more monolithic. The question remains regarding why the cowl's shielding 
performance is at best only moderate. 

We believe these duct diffraction calculations to be correct, if for no other reason because 
the bo undar y-value problem of the rigid cowl as purely "lifting" object is a fairly straightforward 
one (as opposed to when the cowl is lined, to be discussed below). All of the numerical results 
reported here have converged with respect to every computational parameter in the lifting-surface 
theory, and, again, the predicted scattered fields contain the physical features expected of them 
(the specular lobe patterns, etc.). 

One promising hypothesis for the lack of a deep shadow in these rigorous calculations of 
cowl diffraction is the ever presence of a pair of modes whose cut-off and cut-on frequencies 
straddle the frequency of interest, e.g., the BPF. The cut-on member of the pair would naturally 
have the greater impact. At conventional cut-on, based on "infinite-duct thinking," there is a 
severe mismatch between the axial impedance of the fluid inside versus that outside the cowl. 

The fluid within the duct behaves then as an axially incompressible air column with radial 
corrugations corresponding to the acoustic wavenumber divided by the Prandtl-Glauert factor 
p^. I.e., the mode's "z" wavenumber becomes by definition zero. The duct's interior becomes 
an effective resonant space with respect to the two straddling modes. 

Figs. 15a,b, respectively from Meyer 7 and the author’s Ref. 8, support that working 
hypothesis. Fig. 15a shows computed reflection coefficients versus frequency for an inlet duct 
driven from within by the fan stage of an Advanced Ducted Propeller having the same number of 
rotor and stator blades as that of our calculations. The reflections are largest at cut-on for the 
reasons just cited, viz., due to the forced coexistence of an axially hydrodynamic interior and a 
standard wave-bearing exterior. What is interesting, as well as important to the present 
argument, is that the amplitudes of the reflections remain fairly high over wide bands of 
frequencies: the cowl's interior leaks energy to the exterior but nonetheless stays largely locked 
in a semi-resonant state. The resistance supplied by the exterior is not enough to keep that from 


34 



CAMBRIDGE ACOUSTICAL ASSOCIATES, INC. 


happening. 

Fig. 15a is actually a recreation of Fig. 15b in a different guise. The solid curve in this 
second figure plots the formally computed power radiated by a single axial dipole within a rigid 
unflanged duct of finite length. There is no freestream. The ordinate has been normalized by the 
dipole's power output in freefield, so that values above zero represent an amplification of the 
leaking acoustic field due to the interior's quasi-resonant state. Fig. 15b's solid curve remains 
consistently above that zero level even at the highest frequencies "k" in the graph (the abscissa k 
stands for oL/2c in the nomenclature of this report). The author further recalls that the resulting 
amplified directivity patterns (not published in Ref. 8) did not always beam away from the duct s 
nominal shadow to collect in the duct's unprotected sectors, so that much of the added radiation 
in fact went to fill in the shape's expected shadow. This comment would likely apply also to 
Meyer's predictions, perhaps even more so because they address relatively high-index spinning 
modes with naturally low near-axis levels (Meyer does not report directivity patterns in Ref. 7). 

To summarize the dual message of Figs. 15a,b: 15a proves that inlet reflections near cut- 
on are high over bands of frequencies that are surprisingly broad; and Fig. 1 5b concludes that 
even thou gh the semi-resonant fluid interior is reactively isolated from the exterior, the same is 
clearly not true of its resistance. The resulting power output of the partly encased dipole is 
greater than it would have been in freefield, over wide bands of frequencies: The true waveguide 
character of the rigid cowl largely nullifies the shielding predicted canonically for the unwrapped 
shape. Section B.2 will explore whether practical arrangements of a lossy liner might 
sufficiently dampen the interior's quasi-resonances and thereby allow the waveguide system to 
approach the ideal canonical state. 

B.2 Calculations of Cowl Diffraction f or a Lined Cowl 
B.2a Preliminary Remarks 

The liner calculations will address only the BPF (m=l). It was earlier found that for the 
stator blades in question Amiet's 21 chordwise compactness parameter produced a typical value of 
.54. That number is lower than the 7t/4 rough dividing line beyond which noncompactness 
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effects are so pronounced that a leading-edge solution is a good approximation to the exact 
inversion of the Possio kernel and its convolution with the input gust 

For no particularly compelling reason, however, the brief liner study about to be 
discussed does not invoke the chordwise compactness condition even though the radiating stator 
blades meet it from root to tip (i.e., ".54«*/4"). The calculations instead treat all blade loadings 
as chordwise noncompact. The stator's incident field will therefore not match exactly those 
earlier computed via the alternate chordwise-compact code. 

All liner cases will use the same ring-band patch. The cases examined will shift that 
patch axially along the cowl's interior to determine the effect of its position on overall cowl 
shielding. The locally reacting microducts that make up the liner's surface will be tuned to the 
driving frequency, i.e., to the BPF: RPM^ in Eq. A.9 will be the rotor's RPM. The loss factor 
H will be unity. The movable lined ring band will measure one-third of the cowl's length L: 

L]iner = L/3 . 

All predictions will be for the ground reference frame relevant to community noise. 

B.2c Formal Diffraction bv a Fined Cowl; Numerical Results 

Fig. 1 6a shows predictions of incident and total fields for an inlet placement of the liner 
segment (see upper engine sketch). The cowl’s shielding performance is somewhat improved 
over that earlier displayed in Fig. 13a for the perfectly rigid case. Fig. 16b compares the incident 
field to that scattered by the now-partially lined system. This second field differs significantly 
from its rigid-cowl version in Fig. 13d. Fig. 16b’s specular lobes appear split into smaller sub- 
lobes, more so in the forward than in the rear direction. The scattered field remains roughly 

unaffected by the inlet liner for angles near 0 o =9O°. 

Figs. 17a,b repeat Figs. 16a,b but for a central lined segment (for manufacturing and 
maintenance reasons, this case is probably only academically interesting). The cowl's shielding 
performance in Fig. 17a appears to be better yet than in Fig. 16a, except for some degradation 

along the strong forward lobe around 0 O =25°. 

Figs. 1 8a, b are for an outlet liner. This configuration does not appear to perform globally 

as well as the other two, although it reduces noise somewhat along the downstream axis as 
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expected. The outlet liner has not changed much the three forward lobes in Fig. 1 8b's scattered 
field upon comparison to what they were in Fig. 13a. That is again according to physical 
expectations: those lobes are chiefly due to a local bounce-off of stator noise striking once the 
nearby untreated inlet. 

Fig. 19 collects the total fields from Figs. 16a-19a in order to showcase their differences. 
The picture's conclusion is that each of the three liner layouts enhances the cowl s rigid shielding 
performance in a roughly understandable way: the inlet liner dissipates some of the forward fan 
noise and so, reduces somewhat the total field within the main forward lobe at 0 O =25°; and (2) the 
outlet liner similarly cuts down the already-low rearward radiation; (3) the center liner sends 
slightly less acoustic energy both forward and backward than the rigid case. The radial content 
of that energy appears to be such, that while it improves the cowl’s shielding performance for 0 O 
near 90°, it does so in part at the expense of fore and aft radiations. 
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VI CONCLUSIONS 

(1) The project has produced a fairly sophisticated model for the cowl's incident field due 
both to stator/rotor and to rotor/stator interactions. The source analysis has relied on strip 
aerodynamics, although all blades and impinging wakes were otherwise three-dimensional. The 
blades were twisted from root to tip and their spanwise-variable chords were aeroacoustically 
wide, or noncompact. The theory has generated formal expressions for both the nearfield 
pressures and fluid particle velocities along the cowl s diffracting walls. 

(2) Part of the work centered on the development of a novel set of circumferentially 
unwrapped, but axially warped canonical models to predict the ideal performance of a cowl of 
arbitrary generator curve under conditions of perfect absorption for its interior. The remaining 
acoustic rays were made to diffract around the ends of a two-edged generalized lh 3 flg w- 
Those canonical calculations determined that cowl camber and/or taper should be inefficient 
agents for steering rotor/stator noise away from near-axis directions. The numerical predictions 
of the canonical models also supplied important benchmark features against which the physical 
content of the later calculations of rigorous cowl diffraction could be judged. 

(3) The formal predictions of cowl diffraction of the rotor/stator incident fields, 
described under (1) above, led to the following conclusions: (a) an unlined system should 
perform at best moderately, whether its insonifying blades are chordwise-compact or 
noncompact, and both for the blade-passage frequency and for its overtone, (b) A lined cowl 
should perform somewhat better, with typical shieldings in the range of several dB. That higher 
performance came to be roughly consistent here with one's expectations: an inlet liner has its 
greatest effects in the forward direction, etc. 
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leading-edge coordinate comes into play in the noncompact version of the airloads 
theory, (c) Rotor-first, stator-second arrangement typical of the fan stage of most 
modern commercial aircraft engines. 
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. Modeled Gaussian form of the mean wake deficit u for a stator blade. Variable $ ranges 
ftoToto 2a IB S . and the distance indicated responds to the U2-«ormahzed form of 

the distance X in Ref. 12'sFig. 3.18. 
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Figs. - 4a,b (a) Original and geometry of ring-wing chord L, insonified by a three dimensional 
volume velocity source of strength Q o at z = z * in Mach number flow . The 

diffracted field affecting a listener directly overhead of the cowl is determined 
primarily by the points A and B on the inlet and outlet, respectively. 

(b) (left) unwrapped duct to provide the locally flat, two-edge model invoked in later 
figures, indicating the positions of points A and B. The implied observer is on a 
plane normal to the flattened "ring wing". The intersection of that plane and the 
wing is the chord line containing the points A, z - z * , and B; (b)'s right-hand sketch 
shows the wing's chordwise view. 
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Figs. - 4c,d (c) Behavior of the diffraction load Lp and the jump in potential A $ on approach 
of the leading-edge of a semi-infinite rigid wing, (d ) Candel and Amiet's trailing-edge 
diffraction problem for a pressure source of strength " P 0 "• A p vanishes on approach 

of the trailing edge as the square root of the distance to it, i.e., as the leading-edge 
behavior of the potential jump A $ in the upper figure. 
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Fig - 6 Coherent combination of the shadow zone fields of Figs. 5a,b, leading to a new fringe 
in that shadow region. The bottom sketch re-wraps the cowl to recall the relative 
position of source and edges in the original geometry. 
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f; „ 7 o b Su b-canonical calculations to verify Figs. 6's directivity pattern: (a) lobe field in the 

Fl g, . 7a,b Sub^ano^ ^ ^ ^ factors oftwo "diffiactron" sources a distance 

L apart in a freestream M k ; (b) The lobe pattern on the illuminated side of the 

leading-edge solution in Fig. 5a, compared to that of a source near anMmte rigid 
wall. The two patterns match towards the right side of the picture as they should. 
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Fit. 1 la- Aberration of sound: necessary tilting of ray directions for achieving a spatial 
obtoive. The figure shows the cowl's Wet Mid outlet stations, and a typical poult of radiation at 
the cowl’s axW center but removed fiom it radially by an aibitrary amount. Rays shooting 
upstream must do so at an "actual” angle <££ or less if they are to clear the cowl (upper 
equation). Downstream rays that make an angle or greater measured from the 

downstream axis strike the cowl. 
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Fie lib : Effect of ray filling on the cowl's downstream zone of specular reflations (Snell's law 
in a flow). The picture shows a typical origin of radiation, and the "last" ray from that pom o 
experience a single bounce against the cowl's inner wall before proceeding to e ower 
medium. The flow endows the downstream zone of specular reflections with a ray-directional 

content greater than that of a still medium. 
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SSSl The rafflating rotor blades are justifiably chordwise compact. The 
cowl is unlined. 
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Fig 13 a: Incident and total fields for a practical rotor-first configuration: same as Fig. 12 but 
with the fen stage reversed. The cowl is unlined. 




Fig 13b Total pressure fields as perceived on the aircraft’s reference frame (solid), and on the 
ground (dashed). The former is washed downstream by the meanflow. The ^e s oimd field 
beams forward when heard at ground level by an observer stationary relative to the air. 
flight Mach number is .4. The cowl is unlined. 
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Fi g 13c : Comparison of the stator’s incident fields for observers riding with the aircraft (solid), 
and at ground level (dashed). 
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Pip i3d : Diagnostic comparison of incident and cowl-scattered fields in the aircraft's reference 
frame. The specular and nonspecular lobes are interpreted as such. The cowl is unlined. 
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Fios 14a.b : Predictions for the overtone m=2 of blade-passage. The twisted radiating stator 
blades* are now chordwise noncompact at every radial station from root to tip. (a)i shows mcident 
and total fields in the aircraft's frame of reference; and (b) shows them as perceived on the 

ground. The cowl is unlined. 
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Fig, I4c.d : (c) Comparison of the incident fields heard when flying with the aircraft and at 
ground level, (d) Similar plot for the total fields. 
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F?g i4e ; Diagnostic comparison of incident and cowl-scattered fields, with the zones of 
specular reflections identified as such. The cowl is unlined. 
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Figs 1 5a_b: fa') From Meyer 7 ; computed inlet reflection coefficients for three duct roodra. 
(ifFromUrtinez«; the solid curve shows the radiated power of a stngfc ^aldipole within 
rigid duct. The ordinate appears normalized by the source s freefield radia p 
medium is still. 









Figs 16a_b: (a) Predictions of incident (solid) and total (dashed) fields for an inlet liner covering 
one third of the duct’s length. The listener is on the ground, (b) Incident and cowl-scattered 
fields. 
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Fi gs. 17 a r b: (a) Predictions of incident (solid) and total (dashed) fields for a central liner 
covering one third of the duct's length. The listener is on the ground, (b) Incident and cowl 
scattered fields. 






Figs. 18 a_b: (a) Predictions of incident (solid) and total (dashed) fields for an outlet liner 
covering one third of the cowl's length. The listener is on the ground, (b) Incident and cowl- 
scattered fields. 
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Fi g. 19: Comparison of the predicted total fields of the three liner configurations in 
Figs. 16-18. 
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APPENDIX A: Coupled Lifting and Nonlifting Surfac e Theories for the Straight Cowl in its 
Role of Aeroacoustir Diflfractor and Dissioator of Proouls or Noise 
(adapted in part from Ref. 1) 

The two fundamental unknowns in the cowl's diffraction problem are the A p y (z) 
diffraction load, which will be nonzero whether or not the duct is lined, and A w v (z) , which is 
nonzero only over the lined segments along the interior side of the cowl's two-sided wall. The 
first of the two coupled integral equations applies to the r=a + outer surface. That tide is taken 
here to be rigid, and the boundary condition is therefore one of flow tangency for the incident and 
the cowl's self-induced field of normal velocities. The negative of the former appears on the left 
of the first of the two coupled equations. The cancelling field of induced flows generates the right 
tide: 



The second equation comes from the boundary condition along r=a* surface, part or all of 
which may be covered with a dissipative liner: 
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where a is the liner's axially variable admittance normalized by the fluid's characteristic admittance 
1/p c: 
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Eq. A.2 collapses to the statement A w y (z) * 0 if there are no lined segments at all along 
the cowl's interior. Eq. A 1 then turns into a more familiar looking lifting problem, uncomplicated 
by the one-sided dynamic thickness effects of the compliant dissipative liner: 
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Kernels K V M , K v * ,lh , and K^* in Eqs. A 1,2 are 
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Complex contour C is as drawn on Fig. A. 1 . Symbol k stands for o)L/2U„, the reduced 
frequency of the cowl based on the cowl's halfchord L/2, and y is the radial wavenumber earlier 
defined with an "m" subscript which has been omitted here for the sake of simplicity: 
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(A.8) 

The liner material is locally reacting and its impedance is that of a blocked acoustic layer 
of thickness h^ and complex compressional wave speed c^*. (the blocking is provided by the 
cowl's "outer" wall, which is of infinite impedance, cf. Fig. A.2). The dissipation mechanism will 
be a material loss factor r| associated with that compressional wave speed. Normalized by the 
medium's characteristic impedance pc, is 




The second equality has fixed the real part of the liner's compressional wave speed to 
make impedance pass through its first null, i.e., through its first thickness resonance, at a 
value of frequency w corresponding to a propeller RPM of choice, RPMa***. 
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APPENDIX B: T ifting- and Nonlifting-Surface Th eori es fpr a Cowl pf Arbitrary Cambe r 
and Taper 


B.l General 

The next two sections document the generalization of Ref. 1 s lift and dynamic thickness 
analysis for a straight duct, or cowl, to now account for fore/aft camber. The cowl’s cross section 
remains circular at every axial station. Section B.2 addresses the steady incompressible limit of 
the theory that would be relevant for applications at low takeoff speeds when the cowl is at an 
angle of attack. The results of that limiting subtheory for a ring wing in steady incompressible 
flow could be used to generate the duct’s field of nonunifonn velocities over the stage. Section 
B.3 extends the new subtheory to cover a more immediate goal: the diffraction kernels for the 
unsteady subsonic problem of interest here. All spatial variables and wavenumbers will be 
dimensional, unlike the main text. 

B.2 Jh£ Tnrnmnressiblft Steady Form o f the Lifting Kernel for a Cowl of General Fprg/ M 
Camber and Taper 

The incompressible form for general camber of Ref. l's Eq. 3 is 

, , d 2 p[ 

V 3 2 D p , < (r,«|) ) z;fl(s),4> > z(s)) - V^dP* + *T"T ” 

to (B.l) 

- [A/(a(I),$,z(s)).<tf] 6[r ~- ( * - MM*)] *(*-*) ’ 

3«(s) a (s) 

where the indicated directional derivative on the right side generalizes the radial derivative on the 
right side of Ref. l's Eq. 3. The tilde on pwill continue to stand for the frequency-domain version 
of the original time-domain A p. The "e" subscript will similarly follow Ref.l in denoting m the 
influence of the single lifting element at s=$. The V superscript stands again for the lifting 
component of the complete lifting/nonlifting theory. The (-) sign on the right side of Eq. 1 is the 


B.l 



same (-) sign on the third term on the right side of Eq. 3 of Ref. 1; i.e., ft in Eq. 1 is ApdS times 
"-5G/dn(S) H , as explained after Eq. 2 of Ref. 1. Here dS stands for a(5)d<t>d$, the element of cowl 
midsurface area at the source position S. 

With 0(§) now as the angle which the local normal at 5 makes with the horizontal (Fig. 
B.l), the generalized directional derivative at the source position becomes 
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One defines the transform pair 
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as on Eqs. 6a,b of Ref. 1, and applies Eq. B.3a to Eq. B.l. Then 
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where H„ (1) is the zeroth-order Hankel function of the first kind and K„ is its modified form. 
Expressing Eq. B.4 in terms of circumferential modes v gives, via Eq. B.3b, 
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Contour C is as portrayed in Fig. A. 1 with p m = 1 , and upon removal of L/2 used in that figure 
as a normalizing coating [ignore for the moment also the resulting pole at k = - <•> / U^. 

The next step toward the new twice nonplanar lifting-surface theory for the diffracting 
cowl is to relate the pressure p e at a field point on the duct's surface to the normal -to-duct fluid 
velocity at that same point. This process requires the intermediate calculation of a set of pressure 
gradients taken along a line, semi-infinite in z, whose landing point is the axial coordinate z(s) of 
that receiver point. Eq. B.5 has anticipated this requirement and has fixed the field s radial 
variable "a(s) n to correspond to the radial value for the control point, while allowing variable z to 
roam a semi-infinite range -°°<z<z(s). 

The direction of the required gradient of p l t (ci(s), <p, z = free) , corresponds to the 
landing point's surface normal »?(s) [Fig. B. 1, bottom left]. This directional derivative is 
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The second half of the just-cited twofold calculation is the conversion of this pressure 
gradient into a normal velocity w(s). This final operation uses the linearized inviscid momentum 
equation for this velocity component w(a(s),4>,z) of the pertuihed flow along direction #?(s) , 
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Eqs. B.6, 7 and 8b now lead to 
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where the thr ee factors within curly brackets in the wavenumber transform on the right side stand 
for the three similarly bracketed terms in the k integrand of Eq. B.6. 


The "-oo" bottom limit of the z and final integral in Eq. B.9a causes that integral to exist 
only for 


Im [£] > Im 
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One finally obtains that 
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Eq. B. 10 states the transfer function linking the flow velocity along direction 1(s) at 
receiver point a(s),4>,z(s) to the elemental force that induces it at source position a(S),<|>,z(s). The 
passage of contour C over the real-axis pole at k — - o> / U m incorporates Eq. B.9b. This 
existence condition holds regardless of the sign which real life could impose on the imaginary part 
of <o; e.g., in an aeroelastic context the sign of Im(co) might have to be allowed to "float". Eq. 
B.9b would put contour C above - to / U in any case and thereby would generate the shed wake 
expected of the lifting element (Ref. 1, p. 15). 

The last steps in the new development here are. (1) Integrate the right side of Eq. B. 10 
with respect to its source arclength variable § and circumferential running angle <J> over the spatial 
domain -S/2<S<S/2, 0<<|><2rc that corresponds to the cowl's twice-warped lifting surface; (2) 
Identify the circumferential modes of the unknown diffraction loading distribution by means of 
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and finally, (3) Enforce flow tangency over the lifting surface by requiring at every s point that the 
total normal self-induced flow cancel the normal component of the incident field. 

The final lifting-surface equation resulting from those three steps is 
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with kernel K, w (s,$) given by 
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Eqs. B.12a,b are the ©(s)*n/2 generalizations of those derived and solved numerically for 
a straight, uncambered cowl in Ref. 1 (Eqs. 19a,b). One notes that K v M (s,S) is not a 
"displacement kernel" of argument s-5, as it was in Ref. 1 with s=z and §=z. Instead its argument 
is now the more general "s,$" pair. 

The second through fifth lines of Eq. B. 12b appear specialized to segments of the duct's 
surface where the running values of the duct radius a(s) are smaller than the radius a(s) of the 
control point s. For the remaining parts of the surface, these four lines would simply exchange 
themselves for their a(s)<a(5) versions, based on the a(s)<a($) part of Eq. B.6. 

The sixth line in Eq. B. 12b, that shows up left-indented relative to the previous four, is the 
negative of their high-k asymptote. This sixth line applies both to a(s)<a($) and to a(s)>a(§). This 
term has been subtracted and added for two reasons, the first mathematical and the second 
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physical: (1) to improve the convergence of the transform’s integrand at the infinite ends of 
contour C; (2) to display the nature of the kernel's expected singular behavior as the running 
arclength 5 approaches the control point's arclength position s. The negative of this term has been 
inverted analytically and appears as the first term on the right side of Eq. B.12b. This high - k 
term collapses, as it should, to the similar spectral subtraction in Ref. 1 for the special case of a 
straight, untapered and uncambered cowl for 6(s)=n/2 for all s. 

B.3 Corresponding Kernels for Unsteady Subsonic Flow 

The forms of the kernels needed for the diffraction problem for a cowl with a general 
fore/aft generator curve follow from analyses similar to that above. The radial wavenumber y 
appears again without its m subscript. Again for a(s) > the generalization of Appendix A's 
Eq. A.5 for 8 (s) * n /2 is 
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The high- it asymptote of the transform's integrand generates a Cauchy term similar to that which 
appears first on the right side of Eq. A.5. The coefficient of that term is now a function of 0 (r) . 
Returning again to Ref. 1, and defining H v (z - z) to be 
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one now finds that the generalization of w ^ for 0 (s) # it/2 is 
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Again a(s) > u(z) in Eq. B.15. Finally, the generalization of K ^ in Eq. A.7 to 0 * it/2 is 
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Fia. B.l Coordinate definitions and geometric parameters^ f or^a 

generally cambered and tapered cowl. The ^T^ole^'marks 
arclength s, which shows a typical virtual dipole . s marks 

the receiver# or control ring station/ which feels p 

aiong ?he focal normal. Angle 6(a) is measured away fro m the 
downstream horizontal plane. The lower left of the ma: in sketch 
shows the — <*<*<s> trajectory of running points * 
which Eg. B. 8 calculates the pressure gradient in the directi 
5(s) . The insert is from Ref. 1/ where the cowl had neither 

camber nor taper • 
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